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ABSTRACT 



Since the introduction of parameter plane techniques by Siljak 
in 1964, much effort has been devoted to the application of these 
techniques to various control system problems. 

This paper concerns the application of parameter plane techniques 
to the design of active filter systems. 

Constant bandwidth curves on the parameter plane are compared with 
frequency response curves to validate the constant bandwidth curves as 
accurate representations of active filter performance. This in turn 
implies that the constant bandwidth curves can be used in the design of 
active filter systems. 

Also presented is an example of the use of constant zeta and omega 
curves on the parameter plane to manipulate the pole locations of 
individual circuits to achieve a desired overall system configuration 
with the proper frequency scaling. 
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1. Introduction 



The use of operational amplifiers in active filters increases the 
flexibility with which signal filters can be designed. Primary reasons 
for this increased design flexibility are the elimination of the need 
for inductors and of the need for impedance matching. Active filter 
circuits can be used over a wide frequency range, oftentimes with the 
added benefits of less cost, smaller size and less weight. ^2j These 
features suggest that active filters provide a suitable answer to a wide 
variety of filtering problems, such as, aerospace medicine data collec- 
tion and the analysis of vibrational fatigue in aircraft structures. 

This paper is a case study based upon an active filter parallel-T 
filter designed by Mr. William Kerwin of NASA’s Ames Laboratory and is 
primarily concerned with establishing that parameter plane techniques can 
be effectively applied to design of active filter systems. Figure 1 is a 
circuit diagram of the filter and Figure 2 illustrates the observed 
experimental filter performance of that circuit. The experimental per- 
formance will serve as the fundamental reference for the evaluation of 
the parameter plane design techniques used in this paper. 

Two general types of parameter plane curves will be considered: 

(1) constant bandwidth curves and (2) constant zeta and omega curves. 

Hollister developed equations for displaying curves of constant 
bandwidth on the parameter plane and defined a constant bandwidth curve 
as: 1^3 j 

A constant bandwidth curve for G(j6J^) = M is a curve drawn 
upon the parameter plane which specifies the relation between 
the parameters necessary if the transfer function G(s), which 
is a function of the parameters, is to have magnitude M at 
the real frequency 
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Figure 2. Experimental Filter Performance 



Nutting developed computer programs to plot the constant zeta and 
constant omega curves on the parameter plane. ^4j If these parameter 
plane curves are frequency normalized, then frequency scaling must be 
performed when cascading several circuits together to form a filter 
system. A technique for proper frequency scaling to allow the manipula- 
tion of system pole locations will be demonstrated. Proper conversion of 
parameter plane root locations into correct frequency scaled values for 
the individual circuits permits the individual circuit transfer functions 
to be multiplied together to obtain the system transfer function which 
represents the physical system. 

In establishing parameter plane techniques as applicable to active 
filter design, the following approach will be made. 

First, a mathematical model for the filter system will be developed, 
that is, the transfer function for the system will be derived. 

Then, a digital computer program for the CDC 1604 will be written 
to evaluate the frequency response of the mathematical model of the 
system. This frequency response will be compared with the observed 
experimental response. If the performance of the mathematical model 
conforms with the actual experimental response then the mathematical 
model will be considered a valid representation of the physical system. 
Furthermore, each of the individual circuits which together comprise the 
overall system will be considered to be accurately represented by their 
respective individually derived mathematical models, i.e., transfer 
functions . 

The frequency response of the transfer functions for the individual 
circuits will then be considered accurate enough to serve as standards 
for the evaluation of the constant bandwidth curves on the parameter 
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plane. Good correlation between the constant bandwidth curves and the 
corresponding frequency response curves then implies that both of the 
curves accurately describe circuit performance and that the parameter 
plane can, therefore, be used as a design tool. 
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2. Mathematical Modeling of the System 

Evaluation of the effectiveness of parameter plane techniques 
could be accomplished in several ways. One method would be the Hard- 
ware approach. Use parameter plane techniques to determine the values 
of selected parameters which will produce a given performance, then 
actually build the system and adjust the chosen parameters to the speci- 
fied values. Obseirve the experimental performance and compare with the 
desired performance to substantiate the validity of the design technique. 

A second method is simulation. In this case a mathematical model 
of the system is developed and the performance of the system is evaluated 
on either a digital or analog computer. \txen facilities are available 
for either analog or digital simulation these methods have the advantage 
of rapid evaluation and excellent flexibility. 

Even if the first method of hardware analysis is applied, the devel- 
opment of the parameter plane curves requires an accurate mathematical 
model of the system. 

Therefore the first step in design using the parameter plane is to 
develop the mathematical model of the system. 

2.1. Derivation of Transfer Functions Through Signal Flow Graph Analysis 
The mathematical modeling of the filter system shown in Figure 1 will 
be based upon the analysis techniques of the signal flow graph introduced 
by S. J. Mason, j^sj The system will be divided into the three circuits 
indicated in Figure 1. Transfer functions for the individual circuits 
will be derived using signal flow graph analysis. The overall system 
transfer function will be the product of the individual circuit transfer 
functions. This overall system transfer function will be evaluated by a 
computer program to be presented in Section 3. 



2.1.1. Signal Flow Graph Analysis of Circuit 1. 




Figure 3. Circuit Diagram for Circuit 1 
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Figure 4. Signal Flow Graph for Circuit 1 
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Non-touching loops taken two at a time, 
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R 2 C 1 



^.1." 2.1.2 Signal-flow Graph Analysis of Active Parallel-T with Resistive 
Feedback. 

The circuit diagram shown below is an active parallel-T network 
and represents both circuit 2 and circuit 3. 




Figure 5, Circuit Diagram, Active Parallel-T Network 
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Signal Flow Graph, Active Parallel-T Network 
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The loops of the signal-flow graph, Lj; 
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Non-touching loops taken three at a time, 

.2 
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then factor Cp + 1) from 



multiply numerator and denominator by 
numerator and denominator 



J5L_ 

k + 1 
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T(S) 
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where 



p = RCS 

Consideration of load impedance, Z, let 
Z = — I — 
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Substitute into equation (2.2-1) 
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2.2. Frequency Scaling the 'transfer Function 

Whenever several circuits are cascaded it becomes necessary to 
frequency scale the complete network to some reference frequency f^ 
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prior to evaluating the overall frequency response. 

In the RC elliptical function filters under consideration the 
’’notch frequency” is 



f„ = 
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The normalized transfer function then becomes 
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Multiply numerator and denominator times a'^ 
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Frequency scaling the transfer function: 

For example, consider filter section 2 
R = 50 K 
C = 638 pf 
k = 2.0 



= 



n R 2 C 2 



C’ = = 0.436 
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K = 2.762 
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To2(S) 
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Frequency scaling the transfer function. 
Similarly for filter section 3 
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Frequency scaling the transfer function. 
For filter section 1 
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For this section the scaling factor, a, is defined; 
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Frequency normalized transfer function for section 1 
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3. Evaluation of the Mathematical Model 



Does the mathematical model developed in Section 2 accurately 
represent the physical system? Can those transfer functions be used 
in developing the parameter plane curves? 

To answer these questions the digital computer program, Program 
SYS RESP, written in FORTRAN 60, was developed to evaluate system fre- 
quency response when the system consists of several active circuits in 



cascade. 
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PROGRAM GYS RESP 

THI:j PKGGIv'AM MAY BE USED TO EVALUATE THE SYSTEM FREQUENCY 
RESPONSE !M F'OTH MAGNITUDE RATIO AND IN DB FOR A SYSTEM 
CONSIST ir'G OF UP TO FIVE ACTIV^ STAGES WHEN THE NUMERATOR 
POLYNO.-U^'L FOR THE TRANSFER Fmm.CtioN OF EACH STAGE CAN RF 
REPRESENT EO <^Y 



AND the DENOMINATOR POLYNOMIAL RY 

A(l) + A(2)*S + A(3)*S*^*-2 + A(A)*S^<•*3 



THE FOLLOWING DATA CARDS WILL RF USED 
CARD 1 FIRST LINE OF GRAPH TITLE 
CARD 2 SECOND LINE OF GRAPH TITLE 
CARD 3 NSTAGE ,NUMPTS 

CARP 4 V(0)» INITIAL VALUE nr OM^n-A 

CARD 5 XK(1),GAIN COEFF I f I Fm T » ST AC,n 



FORMAT ( 6A8 ) 
FORMAT (6A8 ) 
FORMAT (21 4' 
FORMAT ( F10.4 1 
FORMAT(F10.4) 



CARD 6 3( 1 ,J) »NUMERATOR COF'^F T f I ^NT S , ST AC-F 1, 



IN ASCENDING ORDERS OF S* FORMAT ( 4F10. 4 ) 

CAkD 7 A( 1,J) , DENOMINATOR Cnr FF I C I <^NTS , ST AGE 1, 

IN ASCENDING ORDfpS OF S FORMA T ( 4F 1 0 . 4 ) 

CARD 8,9,10 XK(2), B(2,J) ANT' A(2,J) RESPECTIVELY 

REPEAT THIS SEQUENCE TO ACCOHODATF THE REQUIR'^D NUMBER 



OF STAGE" 



DIMENSION X(S),B(5,5),A(5,5), TTITLE(12), W(400), 

1 BREAL( 5,400) ,BIMAG( 8,400 ) , 3M AG ( 8 , 400 ) , DR(5,400), 

2 AIMAG( 5,400 ) , AMAG ( 5 , 40R ) , XK ( ^ ) , A° S ( 900 ) , ARE AL ( F , 400 ) , 
3DBSYS(RC0 ) 

DBSYS(0)= 0.000 
DBSYS (1) = 0.000 
W(0)= 0.000 
W(l) = C.OOO 
LABEL = 4H 

READ 200, ( ITITLEt I) , 1=1,6) 

READ 200, (ITITLE(I), 1=7,12) 

READ 30, NSTAGE, NUMPTS 
READ 5C, W(0) 

DO 60 J = 1, NUMPTS 
ABS(J) = 1.0 
60 CONTINUE 

DO 13 I = 1, NSTAGE 
READ ’0, XK ( I ) 

READ 100, ( B( I ,K) , K = 1 ,4) 

READ 100, (A( I ,K) , K= 1,4) 

PRINT 130 ,XK ( I ) , ( B ( I ,K ) ,K=1 ,3 ) 
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PRINT 13I',(A( I ,K) ,K=1,4) 

DO 25 J = 1,NU.MPTS 
A(J) = W(J-l) + 0.01 

BR1£AL(I,J) = B(I»1) - B(I,3) *W(J)**2 

BIMAG( I ,J ) = B( I »2 )*W ( J ) 

BMAG(I,J) = SQRTF(0RFAL( I + B I f^AG ( l , j ) **2 ) 

AREAL (I,J) = A(I,1) - A(I,3)* W(J)**2 

A I MAG ( I, J) = A(I,2)* W(J) - A(I,4) W(J)**3 

AMAG (I.J) = SQRTF (AREAL (I»J)**2 + A I MAG ( I , J ) **2 ) 

TEMP = X;<( I )*BMAG( I»J)/AMAG( I ,J) 

ABS(J) = ABS(J)*TEMP 
DB(I»J) -- 20.*LOG10F( TEMP) 

25 DBSYS(J) = DBSYS(J) + DB(I,J) 

13 CONTINUE 

DO 27 J = 50,150,5 
PRINT 310, (DB( I ,J) ,1=1,3) 

PRINT 320, DBSYS(J) 

27 CONTINUE 

CALL DRAW(NUMPTS,W,DeSYS,0,0,LAPEL,lTlTLF,0,10. ,-^,1 ,0,2,8, 
1 12,1, LAST) 

CALL DRAW (NUMPTS,W,ABS, 0,0, label, I TITLE, 0,0, 0,1, 0,2,8, 

1 12,1, LAST) 



30 


FORMAT 


(214) 






50 


FORMAT 


(F10.4) 






IOC 


FORMAT 


(4F10.4) 






130 


FORMAT 


(1X,30HTHE XK 


AND S coefficients ARF 


,/ ,4F10.4) 


150 


FORMAT 


(1X,30HTHE a 


COEFFICIENTS ARF 


,/,4El0.4) 


200 


FORMAT 


(6A8 ) 






31C 


FORMAT 


(5X,18HDB OF 


THE STAGED ,3E18.6) 




320 


FORMAT 

END 

END 


( 5X, 18HDBSYS 


EQUALS ,E18.6) 





JOB 0600 NAKAGAWA PARALLEL-T ACTIVE FILTER 
X = W(J) Y = DB OUT SK = 2.0 

3 400 
.00 
.9316 
1.000 

.4658 1.1553 l.CCO 



1.1967 






2.4491 


0.000 


1.000 


1.0611 


0.112 


1.000 


1.1198 






1.5220 


0.000 


1.000 


.822 


0.493 


1.000 
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The system frequency response plotted by Program SYS RESP when 
the previously derived, frequency scaled circuit transfer functions are 
inserted as data is shown in Figure 7. Since it can be shown that 3.18 
kc corresponds to a normalized omega equal to unity, comparison of this 
normalized system response with the experimental filter performance 
recorded in Figure 1 clearly indicates that the mathematical model is 
valid, that is, the transfer function is correct. In addition these 
results also indicate that the program will accurately simulate the per- 
formance of the physical system. 

Validity of the system mathematical model implies: 

(1) that each of the circuits which together comprise the system 
evaluated above can in turn be represented by its respective transfer 
function. 

(2) that the performance of these individual circuit transfer 
functions can now be used in turn to evaluate the constant bandwidth 
curves on the parameter plane. 



27 



Amplitude Ratio, db 



\ 



I 




Figure 7. Digital Simulation of System Frequency Response 
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4. Active Filter Constant Bandwidth Curves on the Parameter Plane 

The computer program included in Appendix A was written in FORTRAN 
63 by Jean Bow, Programmer, Computer Facility, Naval Postgraduate School 
and modified slightly by the author to display constant bandwidth curves 
on the parameter plane. The program is based upon the constant bandwidth 
equations developed by Hollister and utilizes the transfer function for 
the active parallel-T network derived in Section 2.1.2. 

The constant bandwitdh curves with amplitude contours in decibels 
are shown on the parameter plane plots of Figure 8-1 through Figure 8-18. 
As noted earlier each parameter plane plot of the constant bandwidth 
curves pertains to a fixed value of 

To obtain a sketch of the frequency response of the circuit, select 
the values of the two parameters on the abscissa and ordinate. This 
establishes an operating point on the parameter plane. Enter each plot 
at the selected operating point to determine magnitude, M, at the fixed 
value of represented by that plot. The frequency response curve can 
now be constructed by plotting M as a function of^. 

Is this technique valid, that is, is it accurate enough to be con- 
sidered as a possible active filter design technique? 

The computer program included in Appendix B was developed specific- 
ally to answer these questions. Written in FORTRAN 60, Program 
FREQRESP evaluates the frequency response of the parallel-T active fil- 
ter. 

Figure 9-1 is the graphical output for this program with SK, the 
ratio of the circuit element values, set at 2.0 and C, the normalized 
capacitance, equal to 0.20. The amplifier gain is varied from 1.6 to 
5.4 to produce the family of curves shown, i.e., each curve represents 
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The following information applies to Figures 8-1 
through 8-18: 

X Variable ’’Amplifier Gain, K” 

X Scale 1 inch = 0.5 

Y Variable ’’Normalized Capacitance, C” 

Y Scale 1 inch = 0.1 
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Figure 8-1 



Active Filter Constant Bandwidth Curves on Parameter 
Plane k = 2.0 ^ = 0.10 
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Figure 8-2. Active Filter Constant Bandwidth Curves on Parameter 
Plane k = 2.0 €U « 0.20 
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Figure 8-3. Active Filter Constant Bandwidth Curves on Parameter Plane 
k = 2.0 0 ) = 0.30 
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Figure 8-4. Active Filter Constant Bandwidth Curves on Parameter Plane 
k = 2.0 U) = 0.40 
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Figure 8-5. Active Filter Constant Bandwidth Curves on Parameter Plane 
k = 2.0 OJ = 0.50 
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Figure 8-6. Active Filter Constant Bandwidth Curves on Parameter 
k s 2.0 CO - 0.60 
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Figure 8-7. Active Filter Constant Bandwidth Curves on Parameter Plane 
k = 2.0 CJ = 0.70 
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Figure 8-8. Active Filter Constant Bandwidth Curves on Parameter Plane 
k = 2.0 CJ =: 0.80 
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Figure 8-9. Active Filter Constant Bandwidth Curves on Parameter Plane 
k = 2.0 CU = 0.90 
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Figure 8-10. Active Filter Constant Bandwidth Curves on Parameter Plane 
k = 2.0 Ci) - 0.95 
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Figure 8-11. Active Filter Constant Bandwidth Curves on Parameter Plane 
k = 2.0 = 1.05 
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Figure 8-13. Active Filter Constant Bandwidth Curves on Parameter Plane 
k = 2.0 CJ = 1.20 
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Figure 8-14. Active Filter Constant Bandwidth Curves on Parameter Plane 
k = 2.0 = 1.30 
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Figure 8-15, Active Filter Constant Bandwidth Curves on Parameter Plane 
k = 2.0 CO = 1.40 
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Figure 8-16. Active Filter Constant Bandwidth Curves on Parameter Plane 
k = 2.0 CO = 1.60 
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Figure 8-18. Active Filter Constant Bandwidth Curves on Parameter Plane 
k = 2.0 C = 0.20 
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X Variable 
X Scale 

Y Variable 

Y Scale 



Normalized Radian Frequency 
1” = 0.5 

Magnitude of Output, M (db) 
1" = 5.0 db 
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Figure 9-2. Active Filter Frequency Response Curves 
k = 2.0 C = 0.436 
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X Scale 
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Figure 9-3. Active Filter Frequency Response Curves 
k = 2.0 C - 0.284 



X Variable 
X Scale 

Y Variable 

Y Scale 



Normalized Radian Frequency 
1" = 0.5 

Magnitude of Output, M (db) 
1” = 5.0 db 
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a specific value of amplifier gain. At small values of co and again 
for Cjl) »1.0 we observe the expected result that output increases as 
gain increases. But the distinguishing feature of this family of curves 
is that near the ”notch frequency” of 1.0 a resonance peak is observed at 
some intermediate value of amplifier gain. 

To evaluate the constant bandwidth curves, two operating points were 
selected. The first operating point corresponds to Circuit 2 in Figure 
1. For this circuit 
SK = 2.0 
C = 0.436 
K = 2.762 

The frequency response curves for this circuit are given in Figure 9-2. 

A comparison of constant bandwidth against the frequency response values 
is provided in Table I. 

The second operating point selected corresponds to Circuit 3 of 
Figure 1 for which 
SK = 2.0 
C = 0.284 
K = 2.074 

The frequency response curves for Circuit 3 are displayed in Figure 9-3 
and tabulated in Table II. 

Other operating points could well have been selected, but further 
repetition would do little more than reemphasize the conclusions that 
are arrived at by considering Tables I and II. For it is clearly appar- 
ent that the constant bandwidth curves on the parameter plane are excel- 
lent predictors of circuit performance. 

It is therefore valid to use the constant bandwidth curves for the 
design of such active filters. 
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EVALUATION OF "CONSTANT BANDWIDTH" CURVES 



TABLE I - Circuit 2 (C = 0.436, K = 2.762) 



RADIAN 

frequency 


AMPLITUDE FROM 
CONSTANT BANDWIDTH CURVES 


AMPLITUDE FROM 
FREQUENCY RESPONSE CLTiVES 








0.1 


9.0 db 


9.0 db 


0.2 


9.3 


9.4 


0.3 


10.0 


10.0 


0.4 


11.2 


11.3 


0.5 


13.5 


13.6 


0.6 


19.5 


20.0 


0.7 


18.5 


19.0 


0.8 


6.0 


6,0 


0.9 


-4.4 


-4.4 


0.95 


-12.1 


-12.5 


1.05 


-14.7 


-14.0 


1.10 


-9.8 


-9.0 


1.20 


-5.7 


-5,6 


1.30 


-3.7 


-3.6 


1.40 


-2.5 


-2.5 


1.60 


-1.2 


-1.1 


1.80 


-0.4 


-0.3 
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EVALUATION OF ’^CONSTANT BANDWIDTH” CURVES 



TABLE II - Circuit 3 (C = 0.284, K = 2.074) 



RADIAN 

frequency 


AMPLITUDE FROM 
CONSTANT BANDWIDTH CURVES 


AMPLITUDE FROM 
FREQUENCY RESPONSE CURVES 




tl 


tl 


0.1 


6.2 db 


6.2 db 


0.2 


6.5 


6.5 


0. 3 


6.8 


6.8 


0.4 


7.1 


7.1 


0.5 


7.5 


7.5 


0.6 


7.8 


7.6 


0.7 


6.0 


6.0 


0.8 


1.8 


1.4 


0.9 


-6.1 


-6.4 


0.95 


-13.7 


-14.0 


1.05 


-15.7 


-15.0 


1.11 


-10.6 


-10.6 


1.20 


-6.3 


-6.5 


1.30 


-4. 3 


-4.3 


1.40 


-3.0 


-3.1 


1.60 


-1.8 


-1.8 


1.80 


+0.1 


-1.0 
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5. Active Filter Constant Zeta and Omega Curves on the Parameter Plane. 

The computer programs developed by Nutting were the basis for the 
constant zeta and omega curves displayed on the parameter plane of Fig- 
ure 10. The transfer functions used for this parameter plane plot are 
derived from the active filter sections that have been designated circuit 
#2 and circuit #3. 

By careful selection of the proper values of amplifier gain, K, and 
the normalized input capacitance, C, it should be possible to adjust the 
root locations to any desired value of damping ratio, zeta, and radian 
frequency, omega, displayed on the parameter plane of Figure 10. 

5.1 Pole Locations and the Design of the Low Pass Filter 

The ideal low pass filter which is frequency response optimized 
has a gain vs. frequency characteristic shown in Figure 11. 

Although such filter performance is not physically realizeable, 
there are many filters that provide acceptable approximations to the 
ideal. Two mathematically designed filters, the Butterworth and the 
Chebyshev, will be discussed here. 

The Butterworth filter is said to be optimally flat at zero fre- 
quency, with amplitude related to frequency by 

1 

A I 2n 

where 

|a| is the amplitude 

is the cut-off frequency 

n is the order of the filter transfer function 

This relationship is plotted in Figure 12. 

The frequency response of the Chebyshev filter shown in Figure 13 
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Normalized Input Capacitance, 




56 



i 

i 



approximates the ideal frequency response with a "constant ripple" re- 
sponse in the pass band, confining the error to some specified value 
throughout the pass band. It should be noted that the Chebyshev filter 
places equal emphasis on all frequencies in the pass band, whereas the 
Butterrorth filter emphasizes the near zero frequency behavior. 

To achieve the frequency response shown in Figure 12, the poles of 

the Butterworth filter transfer function displayed on the s-plane are 

uniformly spaced on a semicircle of radius Cd as shown in Figure 14. 

c 

In a similar fashion the frequency response of the Chebyshev 
filter is accomplished by placing the poles of the filter transfer 
function on an ellipse on the s-plane. It has been shown that the proper 
positioning of the poles can be related to the poles of a Butterworth 
filter of the same order by reducing the real component of the respective 
roots by the factor tanh ^ ’ without modifying the imaginary component. 

Thus 



Chebyshev ^^^utterworth^ 



where 



n 



sinh 



-1 1 



n is the order of the filter 

G is related to the amplitude of the ripple 
1 



with 



the minimum amplitude in the 



\ZTTeT 

pass band as shown in Figure 13. 

5.2. Manipulation of Circuit Poles into Chebyshev Configuration 

The poles of the active filter system shewn in Figure 1 will now 
be modified utilizing parameter plane techniques to form a Chebyshev 
configuration of poles. 
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Figure 11. Ideal Low Pass Filter Frequency Response 




Figure 12. Butterworth Filter Frequency Response 




Figure 13. Chebyshev Filter Frequency Response 
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Locus of Butterworth roots 
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Figure 14. Relative Position of Poles for Butterworth 
and Chebyshev Filters (6^^ Order) 
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The purpose of this exercise is two-fold: 

(1) To illustrate the methods required for proper pole relocation 
which must include an adjustment of frequency scaling factors introduced 
through the normalization of the notch frequencies that was one of the 
steps in producing the normalized parameter plane of Figure 10. 

(2) To demonstrate that the constant zeta and omega curves on the 
parameter plane are an effective means of obtaining a desired performance 
from an active filter system. 

For the purposes of this discussion the roots associated with the 
low pass filter circuit #1 will be assumed to be fixed and the positions 
of the remaining poles will be manipulated to obtain the Chebyshev con- 
figuration. Another assumption is that the notch or "zero output" fre- 
quency is high so that the effects of the transfer function zero will be 
ignored, that is, only the poles will be considered. The two assumptions 
are made for the convenience of clarifying the following presentation. 

No serious limitations upon the design of active filter systems are 
introduced by these assumptions. 

Figure 15 is an s-plane plot of the pole configuration of the active 
filter system. The nimibers in parentheses identify the circuit number 
which is associated with that particular pole pair. Thus the poles 
nearest the negative real axis designated (1) are associated with cir- 
cuit and form the basis for the placement of the remaining two sets 
of poles. 

Since the imaginary component of this pole pair is equal to the 
imaginary component of the corresponding Butterworth filter, the poles 
designated (1) establish the radius, of the reference Butterworth 

semicircle. 
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Figure 15. Active Filter System Pole Configuration 
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Co a) = 


CO^ sin 165° = 0.3638 


II 

o 


= 1.406 

sin 165° 



Similarly the real component of the poles designated (1) establish 
the magnitude of the ripple by specifying the value of tanh 



tanh|5 = 


CHohebyshev _ (D 

^^Butterworth 165° 


tanh = 


-0.5776 _ 0.4254 

(-0.9659) (1.406) 



Now that Co ^ and tanh^ have been established, the remaining pole 
positions can be determined. 



^(3) = 


CO^ sin 135° = 0.994 


^l)(2) = 


CO sin 105° = 1.358 

O 


^ ^ Cheby shev 


tanh^Cr(l)Butterworth = tanh^ cos 165° 

-0.578 


^^^^Chebyshev 


0.4254 60 cos 135° = -0.423 
o 


^ ^^^Chebyshev 


0.4254 (O cos 105° = -0.155 
o 



The values of the damping ratio for the respective roots, zeta, 
will be used on the parameter plane of Figure 10 because of their qual- 
ity of maintaining the same value completely independent of any frequency 
scaling or frequency normalization that must now be performed in order 
to take advantage of the information contained on the parameter plane 
plot of Figure 10 which, it must be recalled, has been normalized to a 
*'notch*’ frequency equal to 1.0. 
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To calculate the values for zeta the following relationship must be 



applied: 



Thus : 





^( 1 ) 



Cheb 



\/^ 



.578 

334 + 0.132 



= 0.847 



^( 3 ) 



0.423 



Cheb 



\fo. 



= 0.392 



179 + 0.988 



^( 2 ) 



Cheb 






0.155 

.024 + 1.843 



= 0.113 



In sijmmary, the following values will produce a Chebyshev config- 
urat ion: 



Circuit 


#1 






^(1) 


= 


0. 


00 


CO (1) 


= 


0. 


364 


cr (1) 


= 


-0 


.578 


Circuit 


#3 






^(3) 


= 


0. 


392 


OJO) 


= 


0. 


994 




= 


-0 


.423 


Circuit 


#2 






^(2) 


= 


0. 


113 


07(2) 


= 


1. 


358 


CT(2) 




-0 


.155 



Since the parameter plane of Figure 10 is normalized to a **notch” 
frequency of 1.0, the calculated frequencies listed above must be divided 
by the transfer function scaling factors, a, discussed in Section 2.2. 
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The scale factors for circuit #3 and circuit were 1.234 and 1.566 



Therefore the parameter plane frequencies, Ci) , are: 

P 






0.9941 

1.234 



0.8056 



^Jp(2) 



1.358 

1.566 



0.867 



To obtain these parameter plane values of omega and zeta, enter 
Figure 10 to determine the required values of normalized input capaci- 
tance, C, and amplifier gain, K. The following parameter values will 
result : 

Circuit ^2 

K(3) = 1.605 a(3) = 1.234 

C(3) = 0.100 
Circuit 

K(2) = 2.04 a(2) = 1.566 

C(2) = 0.107 

Refer to Appendix C for algebraic verification that a Chebyshev 
configuration has been formed by the parameter values selected above. 

5.3. Evaluation of Response of System with Pole Configuration Derived 
Through the Parameter Plane. 

What is the effectiveness of the parameter plane in realization 
of the **desired performance’^ which in this case was constant ripple in 
the pass band? Program SYS RESP, presented in Section 3 was used to 
conduct this evaluation. 

The overall system transfer function required as input data to 
Program SYS RESP demands proper frequency scaling of the individual 
circuit transfer functions according to the formula derived in Section 2: 
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To(p) = 



K(p^ ah 



(3C + l)p2 + a (3.0 - 1.5K) + 3C p + 



Circuit #3 



1.605(p^ + 1.522) 



1.30p^ + l.lOp + 1.522 



Circuit ^2 




5 

1.321p + 0.408p + 2.449 



2.04(p^ + 2.449) 



The fixed pole associated with circuit retains its original 
circuit transfer function: 



The results of the Program SYS RESP evaluation of the active filter 
system designed through the application of constant zeta and constant 
omega curves on the parameter plane are shown in Figure 16. These re- 
sults indicate that the design objective of obtaining constant ripple 
performance has been attained. 




.9316 



p2 + 1.1553P + .4658 
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Output Amplitude, 
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6. Conclusions 



It has now been established that the constant bandwidth curves on 
the parameter plane and the constant zeta and omega curves on the para- 
meter plane provide valid information which may be used in the design 
of active filter systems. 

The constant bandwidth curves can provide a very accurate indica- 
tion of filter performance which is directly related to parameter values. 
In addition the constant bandwidth curves provide the filter designer 
with a good qualitative indication of parameter values which should be 
avoided because of extreme sensitivity of filter performance to varia- 
tions in those specific parameter values. Such regions of increased 
sensitivity are easily recognizeable by the increased density of the 
contour lines on the constant bandwidth plots. 

When desired filter performance can be realized with a known pole 
configuration, the constant zeta and omega curves provide a very effec- 
tive means of achieving that desired performance through the proper 
specification of the various parameter values. These curves can also 
provide the user with information pertinent to the tuning of the filter 
when performance is affected by tolerances of the circuit elements or 
the aging of components. 

It should be noted, however, that the frequency normalized para- 
meter plane can yield meaningful information only with proper frequency 
scaling. Considerable care must be exercised to correctly frequency 
scale the transfer functions when several circuits and/or systems are 
cascaded to achieve some desired result. 

Subject to these conditions parameter plane techniques are useful 
tools in the design of active filter systems. The information contained 
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on various parameter plane presentations should lead to much greater 
insight into the problem of correlating parameter values, pole locations 
and filter performance. Used in conjunction with the tremendous flexi- 
bility and freedom that prevails in the relocation of active filter poles, 
these parameter plane techniques should be substantial assets in further 
development of filter design theory. 
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J w 



• »JOB0600F » r'jAjv/K'M'i'i '■> >G»R» 

PROGRAi>' : oD 

C THIS PROGRA.v, WAS DEVELOPED TO -VALUATE TH- FREQUENCY 

C RESPONSE OF THE PARALLr-L-T ACrrv- PILT"'' WITH A tRANS''-'’ 

FUNCTION HAVING A NUVERATCR POLYNOMIAL OF THE FOP‘'' 

< X- (P>^;-2 + 1) 

AND A DLTU-ilNATOR POLYNOMIAL OP THE FORM 

C 

C ( (SK + 1 ) ”-C+l ) -'-P^^*2 + ( ( SK+1 ) * ( 2-.<v ) /5IE+ ( SK+1 ) *C ) *P + 1 

V- 

c where 

C SK IS THE RATIO OF ELEMENT VALUES 

C C IS THL NORMALIZED CAPACITANCE 

C 

C EACH PLOT WILL 3E FOR A FIXED VALUE OF SK AND C, BUT 

C A FAMILY OF GURVF5 WILL RM PLotTfo jq ^:?FpRp5P^^J VARIOUS 

C VALoES Or AMPLIFIER GAIN K. 

f- 

DIr'.KNS ION W ( SOU ) ,LA( AO ) »XK ( AO ) , I T I TLE( 12 ) ,A ( 5 ) »B ( = ) » 
IdREALC 500 ) »aIMAG( 500 ) ,oMAGN ( 50-^ ) , AREAL ( 5CC ) , AIMAG( 500 ) » 

, 2DB( 500) > AMAGNI 500) 

READ 20u , (IT I TLE ( I ) , I =1 »6 ) 

READ 200, (ITITLE(I), 1=7,12) 

LA ( 1 ) = AH<1 .& 

LA(2) = Ah<l,8 
LA(3) = AHKZ.O 
LA (A) = AHK2.2 
LA(5) = AHK2.4 
LA(6) = AHK2.6 
LA(7) = AHK2.8 
LA(8) = AH<,3.0 
LA(9) = AHK3.2 
LA(IO) = AHK3.A 
LA ( 11 ) = AHK3.6 
LA( 12) = AHK3.8 
LA(13) = AHKA.C 

LA( lA) = AHK4.2 
LA ( 15 ) = AHK4.4 

LA ( 16 ) = 4HI44.6 
SK = 6.0 
C = 0,60 
MAXK = 16 
GAIN = 1.4 
W(l) = 0.000 
NUMPTS = 200 
DO 100 I = 1, MAXK 
GAIN = GAIN + 0.2 



77 



XK ( I ) = GAIN 
B( 1 ) = 1.0 
B(3) = 1.0 
All) = 1.0 

A(2) = ( 2.0-XK ( I ) ) *( SK + 1 .0 ) /S< + (SK+1.0)%C 

A(3) = l.u + (SK,+ 1.0)*C 

DO 25 J = 2»2U0 

>\'{J) = W(J~1) + O.jI 

BREAL(J) = B(l) - 0( 3 J)*^^2 

BIHAG(J) = d(2)*W(J) 

BMAGN(J) = S0RTF(3REAL(J)**2 + BIMAG(J)**2) 

AREAL(J) = A(l) - A( 3 )«-W( J)**2 
A I MAG ( J ) = A ( 2 ) *W ( J ) - A ( 4 ) ( J ) **• 3 

AMAGN(J) = SGRTF(AREAL( J)*^;-2 + A r\AG ( J ) **2 ) 

TEMP = Xk ( I )*6MAGN ( J ) /AMAGN ( J ) 

DB(J) = 8.68569^ LOGF(TEMP) 

25 CONTINUE 

IF (I- 1 ) 33,33,35 

33 MODE = 1 
GO TO 55 

35 IF (MAXR - I ) 37,37,36 

36 MODE = 2 
GO TO 65 

37 MODE = 3 
GO TO 55 

55 CALL DRAW(NUV'PT5,W,DB,MODE,0,LA( I ) , ITItl-,.2 ,3., 8, 
10, 2, 0,9, 15,1, LAST) 

100 CONTINUE 
200 FORMAT (6A8) 

END 

END 

JOB 0600 NAKAGAWA PARALLEL-T ACTIVE FILTER 
X=4(J) Y=DB OUT SK=6.0 C=0.60 
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APPENDIX C 



Algebraic Verification of Chebyshev Pole Configuration 

The parameter plane techniques of Section 5 should yield parameter 
values which will result in the desired Chebyshev pole configuration. 

An algebraic evaluation of the characteristic equations for circuits 1, 

2, and 3 will now be conducted to confirm that such a Chebyshev pole 
configuration has, in fact, been attained. 

The parameter plane produced the following frequency scaled transfer 
functions for circuits 1, 2 and 3 respectively: 






0.9316 

+ 1.1553p + 0.4658 



To2^P> 



2.04(p2 + 2.449) 



1.321p^ + 0.408 + 2.449 



T (p) = 1.605(p^ H- 1.522) 

1.30p2 + l.lOOp + 1.522 



The characteristic equations (denominator polynomials) for these 
transfer functions are of the form 

Ap^ + Bp + C = 0 



. B p , £ , 0 

A A 



Thus the natural frequency, 60 , for the equation is: 



OJ 

n V A 



since 
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2k CO^ 



B CO 
^ n 

2C 



and the imaginary component of the pole becomes 






Circuit 1 



0) (1) = 0.4658 = 0.6825 

n 



^( 1 ) = 



(1.155)(0.6825) 

2(0.4658) 



0.847 



COO.) = (0.6825) (0.5331) = 0.364 



Circuit 2 



0 ) ( 2 ) 

n 



2.449 

1.321 



= 1.362 






. (0.*08) (1.362) , 

2(2.449) 



(^(2) = (-1% 362) (0.9936) = 1.354 



Circuit 3 



<y„(3) 



1.522 
1. 30 



= 1.082 






, (1.100)(1.082) ^ 
2(1.522) 



6J(3) = (1.082)(0.9204) = 0.996 



Comparison of these calculated values of zeta and omega with the 
desired Chebyshev values of Section 5.2 substantiates the effectiveness 
of the constant zeta and omega curves on the parameter plane for the de- 
sign task of obtaining a desired active filter system pole configuration. 
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IS. ABSTRACT 



Since the introduction of parameter plane techniques by Siljak in 
1964, much effort has been devoted to the application of these techniques 
to various control system problems. 



This paper concerns the application of parameter plane techniques to 
the design of active filter systems. 

Constant bandwidth curves on the parameter plane are compared with 
frequency response curves to validate the constant bandwidth curves as 
accurate representations of active filter performance. This in turn 
implies that the constant bandwidth curves can be used in the design of 
active filter systems. 
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Also presented is an example of the use of constant zeta and omega 
curves on the parameter plane to manipulate the pole locations of 
individual circuits to achieve a desired overall system configuration 
with the proper frequency scaling. 
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